Integral and derivative dispersion relations in the 
analysis of the data on pp and pp forward scattering. 
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Abstract 



Integral and derivative dispersion relations (DR) are considered for the pp and pp forward 
scattering amplitudes. A new representation for the derivative DR, valid at lower energies 
than the standard one, is obtained. The data on the total cross sections of pp(pp) interaction 
as well as those on the ratio of the real part to the imaginary part of the forward amplitude 
I/-) ! are analyzed within various forms of the DR and high-energy Regge models. It is shown 

that three models for the pomeron, simple pole with intercept larger than one, triple pole 
pomeron and double pole pomeron (both with intercept equal to one) lead to practically 
equivalent descriptions of the data at y/s > 5 GeV. It is also shown that the low-energy part 
of the dispersion integral (from the two-proton threshold up to y/s = 5 GeV) allows one to 
reproduce well the data on p at lower energies without additional free parameters. 
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1 Introduction 



The energy dependence of the hadronic total cross sections as well as that of the parameters 
p = dteA(s, 0)/^smA(s, 0) - the ratios of the real to the imaginary part of the forward 
scattering amplitudes - was widely discussed quite a long time ago (see [H El H] and 
references therein). However, in spite of recent detailed investigations on the subject, the 
theoretical situation remains somewhat undecided, mainly because of the p parameter. 

In the papers [5], all available data on a to t{s) and p(s) for hadron-hadron, photon-hadron 
and photon-photon interactions were considered. Many analytical models for the forward 
scattering amplitudes were fitted and compared. The ratio p was calculated in explicit form, 
from the imaginary part parametrised explicitly by contributions from the pomeron and 
secondary reggeons. The values of the free parameters were determined from the fit to the 
data at s > s min , where y/s min = 5 GeV. Omitting all details, we note here the main two 
conclusions. The best description of the data (with the minimal x 2 /dof, where dof is the 
number of degrees of freedom) is obtained for the model with a to t rising as log 2 s. The model 
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with <Ttot(s) oc s e ,e > was excluded from the list of the best models (in accordance with 
COMPETE criteria, see details in 0). 

Analysis of these results shows that they are due to a poor description of p data at low 
energy. On the other hand, there are a few questions concerning the explicit Regge-type 
models usually used for analysis and description of the data. How low in energy can the 
Regge parametrisations be extended, as they are written as functions of the asymptotic 
variable s rather than the "Regge" variable cos 6t (= E/m in the laboratory system for 
identical colliding particles)? At which energies can the "asymptotic" normalization 

(T t ot(s) = -QmA(s, 0) (1) 
s 

instead of the standard one 

ff " (s) = i SmA(s,0) (2) 
be used? And last, how much do the analytic expressions for p based on the derivative 
dispersion relations deviate from those calculated in the integral form? 

In this paper, we try to answer these questions considering three pomeron models (simple 
pole in the complex- j plane with intercept above unity, double pole and triple pole with 
intercepts equal to unity) for pp and pp interactions at t/s > 5 GeV. 



2 Integral and derivative dispersion relations. 

The amplitude A(s, t, u) is an analytic function of its variables. Consequently it must satisfy 
the integral dispersion relation (IDR) which can be derived from the Cauchy theorem on 
analytic functions. Generally, for proton-proton and proton-antiproton amplitudes, two 
subtractions must be made. However, assuming, in accordance with many analyses, that the 
odderon does not contribute asymptotically, one can show that the dispersion relations for 
pp and pp amplitudes can be reduced to those with one subtraction constant [l]: 

oo 

b e r 

P±(T± 2m p p up J 

where m p is the proton mass, E and p are the energy and momentum of the proton in the 
laboratory system, and B is a subtraction constant, usually determined from the fit to the 
data. The indices +(— ) stand respectively for the pp and (pp) amplitudes. The standard 
normalization (J2J) is chosen in Eq.(j3J). 

In the above expression, the pole contributions and the part of the integral over the 
unphysical cut from the two-pion to the two-proton threshold are omitted because they 
are < 1% (see, e.g. |6J in the region of interest (-y/i > 5 GeV). In other words, we shall 
investigate here only the contribution of the physical region to the dispersion integral. 

The derivative dispersion relations (DDR) were obtained 0] separately for crossing- 
even and crossing-odd amplitudes 

/±(s,0)=A h (s,0)±A_(s J 0). (4) 



E'VE' - E) E'(E' + E) 



p'dE' 



(3) 
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They are very useful in a practice due to their simple analytical form at high energies, 
E > m v : 



$tef + (E,0) « (E/m p ) a tan 



7T / „ d 



Smf+(E,0)/(E/m p 



(5) 



However it is important to estimate the corrections to these asymptotic relations (0) if one 
is to use them at finite s. 

The starting point is the dispersion integral relation for an even amplitude / + (with poles 
and subtraction constant omitted — they will be re-introduced at the end of calculations — ). 
The amplitude is normalized so that 2m p pi a b(J^ = ^smf + (E, 0). Let us consider the relation 
at t = 0: 

oo 

' ' -$mf + (E',0) (6) 



2E 2 

&ef + (E,0) = P 



7T 



E'{E' 2 - E 2 ) 



and represent it in the form 



IE 2 f dE'(E'/m p ) a 
J E'(E' 2 - E 2 ) 



g(E') 



where g(E') = Qmf + (E', 0)/(E'/m p ) a and a is an arbitrary constant (it is not an intercept!). 

Integrating by parts (and using the notations F(E') = (E' /m p ) a ~ 1 g(E'), F'(E) 
all functions), we obtain, assuming that 3m/ + (m p ,0) = 0: 



dF 
dE 



for 
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Now let us transform the integrand functions to the forms 



F'(E') 
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(E'/mp 
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e «'-0/a + e -«'-0/2 

U«'-0/2 _ P -«'-0/2| 



ln|coth-(£'-£)i, 



where £' = \n(E'/m p ) and £ = \n.(E/m p ). 

Taking into account the above expression, one can write the integral © in the form 



Mef+(E,0) 



In 
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(7) 



Rewriting the integrand logarithm (at x = £' — £ 7^ 0) in (JJ|) 
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coth -x 
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and expanding other factors in the integrand of (J2J) in powers of £' — £ (it does not depend 
on £, this is just a trick allowing explicitly to perform the integration and to present the 
result as a series) 4 



where _D ? = J| and £ 7^ 0, 



one can write 



gv _ gv ^gv 

n! 

n=0 



3te/ + (i?, 0) = ^- e C-^ £ -A- £ jr ■ m P, fc, n) ■ ^(0 

7rm p fr^ + 1 7— ^ ^— ^ rein! 



where 



00 



p=0 J fc=0 n=0 







[r(A; + n + 1) + (-l) fc +" 7 (A; + n + 1, 0(2p + 1)] 

and 7(a,x) is the incomplete gamma function. We should note that one can prove that the 
asymptotic expression for DDR as well as the corrections to them do not depend on the 
auxiliary parameter a (thus it is unreasonable to determine the parameter a from a fit to 
data). So we may put, for convenience, a — 1. The first term in (jHj) gives the asymptotic 
form © El H written at a = I. 

After some transformations we obtain for the corrections the series expansion in powers 
of m P l E 

v p=0 

where 

It should be noted that despite an apparent dependence on £ = ln(E/m p ), the above ex- 
pression for C+(p) in fact does not depend on E. This can be proven using the properties 
of exp(-££> c ). 

Then, collecting all terms and adding a subtraction constant B + , we obtain the final 
expression for $lef + (E, 0) 



9fce/+ (£,()) = 5+ + £tan 



1 1 ,-,(-,■>, 00 

7r „ a 



2 di' 



Sm/ + (£,0) 2^ C + (p) /m p yp 

£ ir^2p+l\E ) ' lj 



p=0 



4 Notc that it is not allowed if £ — > 0. 
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A similar expression is obtained for the crossing-odd part of the amplitudes. If a contribution 
of asymptotic odderon is absent then 



9fte/_(£,0) = -Scot 



2dE 



Sm/_ (E, 0) 2 ^ c_(p) /nip \ ^ 



E 



p=0 



where 



CL(p) = |_l_ Sm/ : (jE ,o)} 



3 Phenomenology. 

Our aim is to compare the fits of three pomeron models with p(s) calculated by two methods. 
The first one is the integral dispersion relation, and the second one is the asymptotic form 
of the derivative dispersion relation with a subtraction constant. We would like evaluate the 
importance of the normalization on the quality of the descriptions and on the values of the 
model parameters. 

There are a few difficulties in an analysis of the available experimental data when the 
integral dispersion relations are applied. 

• The imaginary part of the amplitude under consideration must be known in the whole 
kinematic region for a given process, starting from threshold. 

• It is evident that the known parameterization of the cross sections, valid at s > s m j n , 
cannot be used for the part of dispersion integral from the threshold s t h to s min . 

• The available experimental data on the parameters p for all processes, pp, pp, n^p and 
K ± p are of poor quality compared to those for total cross sections. 

In this paper, we apply as a first step the IDR and DDR only to pp and pp data. We 
have fitted the high-energy models to the data at y/s > 5 GeV. 



3.1 Low-energy data. 

Low-energy total cross sections for (143 points) and pp (220 points) interactions at yfs < 5 
GeV were parameterized as functions of the proton momentum in the laboratory system by 
different expressions for different intervals of momentum. 

The values of the free parameters are determined from a fit to the data |8 at s < s m i n 
with the only constraint that the cross sections calculated from a fit to the data at s < s m i n 
must coincide at s = s m i n with those given by the fit to a specific high-energy model, valid 
for s > s min (all data on a and p are taken from [8 ). 

Thus we perform an overall fit in three steps. 

1. The chosen model for high-energy cross-sections is fitted to the data on the cross 
sections only (without p data) at s > s min . 
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2. The obtained "high-energy" parameters are fixed. The "low-energy" parameters are 
determined from the fit at s < s m i n , but with 0"S(s m m) given by the first step. 

3. The subtraction constant B + is determined from the fit at s > s m i n with all other 
parameters kept fixed. 

Then, without fitting, we calculate the ratios p pp and p pp at all energies above the physical 
threshold. The results of such procedure (for the cross sections) are given in Fig.l and Fig.2, 
and will be detailed in a forthcoming paper. 



3.2 High energy. Pomeron models. 

As an example, we consider three models leading to different asymptotic behaviors for the 
total cross sections. For each model, we investigate how the integral and derivative dispersion 
relations work. We start from the explicit parameterization of the total pp and pp cross- 
sections, then, to find the ratios of the real to imaginary parts, we apply the IDR making 
use of the second method for a calculation of the low-energy part of the dispersion integral. 
Then we compare results for ratios calculated through the DDR. 

All the above-mentioned models include the contributions of pomeron, / and to reggeons 
(we consider these reggeons as effective ones because it is not reasonable to add other sec- 
ondary reggeons provided only the pp and pp data are fitted.) 

a% = V(E)+R f (E)±R„(E), (12) 
where Rf = R+, R^ = R~ and 

/ E \<*±(0)-l 

The parameter A can play a role only for the pomeron term in the triple pole model (see 
below), in the simple pole and dipole models A = 1 but we keep it to have a common form 
and notation for the three models. 

When the imaginary part of amplitude is integrated in IDR we consider (for comparison) 
two kinds of normalization: the standard one defined in Eq. (J2J) and the asymptotic one 
given by Eq. For the latter case, in the expressions for cross-sections, the replacement 
E/Xrrip — > s/si is made. 

Besides, we compare our results with the models which are written as functions of — is 
and with the asymptotic normalization (0). We denote such fits as "—is" fits. 

3.2.1 Simple pole pomeron model (SP). 

In this model, the intercept of the pomeron trajectory is larger than unity j7j 



"P 
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We present the results of the fits using derivative dispersion relations (with a subtraction 
constant) and of standard fits with amplitudes defined in accordance with "—is" rule: 

-Af p (s, 0) = ig(-is/ Si r^)-i + ig + (-is/s i r+M- 1 ± g^-is/s.r-^' 1 . (15) 
s 

3.2.2 Dipole pomeron model (DP). 

The pomeron in this model is a double pole in the complex angular momentum plane with 
intercept a-p(O) = 1. 

V(E)= gi + g 2 \n(E/\m p ). (16) 

For the " —is" fit we use 

-A%{s, 0) = i 9l + ig 2 \n(-is/ Sl ) + ig + (-is/s 1 ) a +W- 1 ± g^-is/s^-^' 1 . (17) 
s 

3.2.3 Tripole pomeron model (TP) 

The dominant term at high energy in this model is the contribution of the hardest j-plane 
Regge singularity allowed by unitarity: the triple pole at t — and j = 1. The form of this 
contribution is 

V{E)= gi + g 2 \n 2 {E/\m p ). (18) 

The prescription " —is" gives 

-A%{ S) 0) = i 9l + ig 2 \n 2 (-ts/ Sl ) + ig+(-is/ 's^^' 1 ± g^-ts/ s,)*-^ 1 . (19) 
s 

4 Results, discussion and conclusions. 

Omitting many details of the fits, as well as the numerical values of the parameters, and 
postponing them for a forthcoming paper, we concentrate now on the main results and 
conclusions. The quality of the fits is presented in the Table, where x 2 '/N p (N p is the 
number of experimental points). In Figs. 3 and 4, we show the curves only for the simple- 
pole pomeron model. For other models the curves are practically the same in the region 
where data are available. 

Table 1: The values of x 2 obtained in the various pomeron models and through the different 
methods for the calculation of the ratio p. 





Standard normalization 


High-energy normalization 


IDR 


DDR 


IDR 


" -is",B = 


Simple Pole 


1.0462 


1.0646 


1.1120 


1.1209 


Double Pole 


1.0532 


1.0454 


1.1319 


1.0793 


Triple Pole 


1.0438 


1.0456 


1.1114 


1.1153 
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As one can see from the Table, all models give good descriptions of the data on a to t and p. 
Evidently, the fit with the integral dispersion relations and with the standard normalization 
is preferable. While the data on a are described with x 2 /d°f ~ 0.91, the data on p are 
described less well, with a x 1 /N p « 1.5. However we think that this occurs because of the 
bad quality of the p data. 
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Figure 1: Cross sections in the SP pomeron model at all energies (left) and at at low energies 
(right). The solid lines are the cross-sections tuned to SP model as described in the text. The 
dashed lines are the continuation of the high-energy SP parameterization to lower energies. 



We would like to insist on the fact that the values of p calculated using DDR deviate 
from those calculated with IDR even at y/s < 7 — 8 GeV (see Figs. 3, 4). It means that 
in order to have more correct values of the p at such energies, one must use the IDR rather 
than explicit analytical expressions from the DDR. 

The neglect of the subtraction constant, together with the use of asymptotic formulae in 
a non-asymptotic domain, may be the source of the conclusion of [Sj excluding the simple- 
pole model from the list of the best models. Inclusion of these (non-asymptotic) terms 
improves the description of p considerably, and may lead to different conclusions regarding 
the simple-pole model. 

However, in order to have these final conclusions, one will have to make a complete 
(a la COMPETE) analysis of the all data, including cross sections and p for np and Kp 
interactions. 
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Figure 2: The ratios p in the SP model. Left: the integral and derivative dispersion relation 
fit with the standard normalization. Right: the IDR fit with asymptotic normalization and 
"—is" -fit (see text for details). 
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